A Sterile Neutrino Search with Kaon Decay-at-rest 
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Monoenergetic muon neutrinos (235.5 MeV) from positive kaon decay-at-rest are considered as 
a source for an electron neutrino appearance search. In combination with a liquid argon time 
projection chamber based detector, such a source could provide discovery-level sensitivity to the 
neutrino oscillation parameter space (Am^ ~ 1 eV^) indicative of a sterile neutrino. Current and 
future intense >3 GeV kinetic energy proton facilities around the world can be employed for this 
experimental concept. 
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I. INTRODUCTION 

A set of experiments and analyses hints towards the 
existence of at least one sterile neutrino. Most notably, 
the LSND \\^ collaboration has reported a 3.8cr excess 
of Ve events from a 17^ decay-at-rest source. The mea- 
surement is indicative of neutrino oscillations around 
lS.m? ^ 1 eV , inconsistent with the solar/reactor [2l-|4] 
and atmospheric /accelerator [5HS| oscillation results and 
the three neutrino picture in general. A "sterile" neu- 
trino, which does not couple to the Z-boson [3], has been 
introduced as a possible explanation for the anomalous 
oscillation signal. 

The LSND result is not alone. The MiniBooNE ex- 
periment has recently probed the Am^ ~ 1 eV pa- 
rameter space with neutrinos and anti-neutrinos. The 
accelerator-based MiniBooNE experiment searched for 
v^ to Vf, oscillations and excluded the two neutrino os- 
cillation hypothesis in the LSND allowed region at the 
98% confidence level jTOl. However, MiniBooNE has also 
probed the parameter space with anti-neutrinos and finds 
Vf, appearance consistent with two neutrino oscillation 
in the Am^ G 0.1 — 1.0 eV^ range. The excess in the 
475 < E,y < 1250 MeV energy range has a 0.5% proba- 
bility of being due to background and is consistent with 
LSND dl]. 

The "reactor anti-neutrino anomaly" is another hint to- 
wards the possible existence of a sterile neutrino \1'2\ [T^ . 
The predicted reactor Vf, spectra, applicable to almost 
all reactor-based neutrino experiments, have been re- 
examined and the ratio of the observed rate to the new 
predicted rate is found to be below unity at the 98.6% 
confidence level. The apparent disappearance of MeV- 
scale reactor anti-neutrinos is consistent with LSND and 
an oscillation signal near Am^ ~ 1 eV^, although the 
result is for Ve disappearance (I'e — ^ V^) rather than the 
related 17g — > 17^ appearance process. 

Models with both one sterile neutrino (3 active + 
1 sterile) and two sterile neutrinos (3+2) have been in- 
voked to explain the anomalies |14H17| . Despite the 
aforementioned results, however, the existence of the 
sterile neutrino(s) is still an open question. A definitive 
measurement of the parameter space favored by LSND, 
MiniBooNE anti-neutrino, and the reactor anomaly is 
required |18j . An experimental idea to probe neutrino 
oscillations near Am? ~ 1 eV^ using a kaon decay-at- 



rest source of v^ along with a large liquid argon time 
projection chamber (LArTPC) detector is presented. 



II. EXPERIMENTAL CONCEPT 



kaons at rest produce a monoenergetic 
^ in 63.6% of decays. A detector placed 



Charged 
235.5 MeV ly, 

0(100 m) from a stopped kaon source can search for 
LSND-like v^ appearance events from v^ (originating 
with the decay-at-rest K^ -^ l-^-^'^fi) oscillation in a nar- 
row reconstructed energy window around the expected 
monoenergetic signal. 

The signature of a i/g appearance event is the 
charged current, quasi-elastic interaction v^n — > e^p at 
235.5 MeV. These signal events compete with the anal- 
ogous v^ interaction v^n — > fi^p as well as background 
charged current interactions of i^e from the three body 
decays K+ -^ n^e+Ve (BR=5.1%) and K^ -^ ir^e^vl 
(BR=40.6%). Although such events can nominally be 
considered a background for the v,, appearance search, 
they can be used for flux and cross section measurements 
beneficial to the oscillation analysis. 

The Booster Neutrino Beamline (BNB) at Fermilab 
is a pulsed source of 8 GeV kinetic energy protons. The 
beamline nominally services the MiniBooNE [TTT and Sci- 
BooNE [20, detectors and will provide protons for Mi- 
croBooNE |2T] beginning in 2014. These experiments 
rely on pion decay-in-flight neutrino production for cross 
section measurements and oscillation sensitivity. The 
BNB is currently providing approximately 3x10^° pro- 
tons on target (POT) per year. The beam delivers about 
4x10^^ POT at 3-5 Hz with a total spill time window of 
1.6 /^s. 

Kaons are produced via a large copper target block 
placed directly in front of the primary proton beam. 
Kaons and other unstable charged particles created in the 
primary proton-copper and secondary interactions come 
to rest quickly and subsequently decay. Copper has been 
chosen due to its high thermal conductivity as well as its 
low radiation length (ATq = 1.43 cm), which significantly 
reduces intrinsic Vf. background from K^ and K^ decay- 
in-fiight. There is no magnetic focusing horn required. 

A 2 kton LArTPC is considered as the neutrino detec- 
tor for this study. The detector is envisioned 160 m from 
the target in a direction opposite or nearly opposite the 



primary beam. This orientation reduces decay-in-flight 
induced I'e background. The basehne has been chosen 
in consideration of sensitivity to the LSND allowed re- 
gion as well as event rate given the 1/r^ dependence 
of the neutrino flux. Note that a baseline of ^240 m 
from the neutrino creation point corresponds to oscil- 
lation maximum for 235.5 MeV neutrinos at the LSND 
best-flt Am^ = 1.2 eV^. A schematic of the experimental 
layout is shown in Fig. [Tl 
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FIG. 1. A schematic of the experimental design. The proton- 
on-target induced kaon/neutrino production location is seen 
on the right and the LArTPC-based neutrino detector is seen 
on the left. The most relevant kaon decay modes and domi- 
nant neutrino interactions are shown as well. 

The BNB can be considered an example beamline for 
the experimental concept described here. There are a 
number of other facilities around the world that can be 
used for such an experiment. The main requirements are 
that the proton beam have >3 GeV kinetic energy and 
that it be high power. Both requirements translate into 
the single statement that an intense source of stopped 
kaons is needed for this design. An exposure of 2x10^^ 
protons on target is used for the oscillation sensitivity 
estimates here. Note that the total 8 GeV proton flux 
after the planned proton improvement plan at Fermilab 
will exceed 2x10"'^^ protons per hour [22^. The detec- 
tor mass and beam exposure have been chosen in order 
to provide 5ct sensitivity to most of the LSND allowed 
region. However, the numbers/parameters employed in 
this study are scalable to other beamlines, detectors, con- 
figurations, and experimental assumptions in general. 



III. FLUX, CROSS SECTION, AND NUCLEAR 
EFFECTS 



The 235.5 MeV v^^ flux and all-energy i/g background is 
determined with the MARS [23_ and Geant4 [24 (Quark- 
gluon string precompound with Bertini cascade model) 
simulation packages. Protons (Tp =8 GeV) are directed 
onto a large pure copper target in the positive z direction. 
Backwards going [cos(0z) < —0.75] neutrinos are consid- 
ered as signal candidates for the purposes of event tabula- 



tion. The loose angular restriction is imposed in order to 
allow some detector placement flexibility-although a "de- 
tector opposite beam direction" policy is optimal. Tar- 
get cooling and a detailed design of the target geometry 
have not been throughly examined in this study. The 
monoenergetic i^^/proton yield is 0.049 and 0.038, ac- 
cording to MARS and Geant4, respectively. The follow- 
ing study employs the Geant4 yield; the MARS result 
can be considered a reference value. The background 
flux prediction, alongside the monoenergetic v^ flux, is 
shown in Fig. [2] It is useful to point out the similarity 
of this experimental idea to that of a neutrinoless double 
beta decay search: an excess of events is sought near the 
endpoint of a well predicted and measured background 
distribution. 
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FIG. 2. The electron neutrino flux expected at 160 m. The 
dominant contribution is from charged kaon decay-at-rest. 
The monoenergetic v^, flux (3.5xl0" i//(3E20 POT)/m^) lo- 
cation is also shown. 



A proper treatment of the nuclear physics affecting 
the initial neutrino interaction as well as the interaction 
products visible to the detector is vital for understand- 
ing the attainable neutrino oscillation sensitivity in this 
design. This experiment depends on the efficient iden- 
tification and precise calorimetric reconstruction of Ve- 
induced charged particle tracks and the nuclear physics 
relevant to neutrino interactions plays a significant role in 
determining how well such events can be reconstructed. 

Neutrino interactions on an argon target are simulated 
using the NuWro event generator |25]. NuWro has been 
chosen due to its use of argon spectral functions in de- 
scribing the momentum and binding energy distributions 
of nucleons within the argon nucleus |26] . A spectral 
function implementation is known to be a more accu- 
rate description of the (argon) nucleus than one based 
on the relativistic Fermi gas model |57|. However, the 
related impulse approximation, in which (1) the neutrino 
interacts with a single nucleon and (2) the products are 
propagated through the nucleus, does not necessarily de- 
scribe the interaction properly, especially at low energy 



transfer [55] . The approximation is a general problem in 
the simulation of neutrino events; the spectral function 
model offers the best chance of the accurate simulation 
of neutrino-nucleus interactions at these energies. 

NuWro also simulates intra-nuclear interactions in 
which the hadron(s) produced in the neutrino interac- 
tion are allowed to reinteract with nucleons in the nu- 
cleus. Such reinteractions can result in multiple ejected 
nucleons/hadrons and ultimately impact neutrino en- 
ergy reconstruction. A description of NuWro's intra- 
nuclear simulation can be found at Ref. |2S]. The GE- 
NIE neutrino event generator [28] has been used as a 
cross check of the nuclear model and the neutrino event 
simulation in general. The v^ and v^ charged current 
cross sections at 235.5 MeV employed in NuWro are 
1.3 X 10""'^ m^/neutron and 1.9 x 10"'*'^ m^/neutron, re- 
spectively. The Ve cross section in GENIE is about 25% 
higher than the NuWro prediction. This is not unex- 
pected, as the spectral function implementation generally 
reduces the (differential) cross section, most noticeably 
at low energy transfer PF]. Note that the multi-nucleon 
component of quasi-elastic-like neutrino interactions at 
these energies |29|[5U] . not simulated here, can act to en- 
hance the cross section by more than 10%. The NuWro 
cross sections across the relevant energy range are shown 
in Fig. d 
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FIG. 3. Electron and niuon neutrino cross sections from the 
NuWro neutrino event generator. The cross section is re- 
ported "per target nucleon (neutron)". 



IV. DETECTOR RESOLUTION 

The nominally distinctive monoenergetic v^ oscillation 
signal is obscured by nuclear physics and detector res- 
olution effects, both working to reduce the (available) 
precision in reconstructed neutrino energy. This v^ ap- 
pearance search requires the precise calorimetric recon- 
struction of the electron-induced electromagnetic shower 
and the associated outgoing proton or protons, recall- 



Proton target 


Copper 


^"^^roLr^^yieldatTp-SGeV 


0.038 


Exposure 


2x10^^ protons on target 


Baseline 


160 m 


Neutrino target 


"'"Ar (22 neutrons) 


Neutrino target mass 


2 kton 


Detection efRciency 


100% 


Ve. o-cc at 235.5 MeV 


1.9 X 10-*^ m^ /neutron 


ffi o-cc at 235.5 MeV 


1.3 X 10"'*^ m^ /neutron 


^P for e~ reconstruction 


0.33/V-B(MeV)-(- 0.012 


^^ for proton reconstruction 


0.10 


Background syst. uncertainty 


25% 



TABLE I. The most relevant parameters employed in this 
study. 



ing that multiple nucleons can be emitted due to intra- 
nuclear effects. The ICARUS collaboration estimates 
that the energy resolution for electromagnetic showers in 
a LArTPC [31 from 50-5000 MeV can be parameterized 
with ^=0.33/v/^(MeV)+0.012. The parameterization 
is used here along with an all-outgoing-proton kinetic 
energy resolution (^) of 10% for Gaussian smearing 
the simulated signal and background outgoing particles' 
energy, imitating reconstruction resolution. The proton 
resolution estimate is loosely based on an ICARUS ex- 
periment exposed to a neutrino beam and their mea- 
sured value of ^ = 3.3% for 50 MeV kinetic energy 
protons [32]. Neutrons are considered undetectable, lost 
energy and few-MeV-scale nuclear de-excitation gammas 
are not simulated in this study. A summary of the 
most relevant experimental assumptions and parameters 
is shown in Table |TJ 

Asymmetric smearing effects due to nuclear physics 
and missing energy from outgoing neutrons created via 
intra-nuclear scattering renders the mean reconstructed 
neutrino energy {E[^^ = ^electron + YTi ^i, proton, where n 
is the number of protons) lower than the true neutrino 
energy (-Etiuc)- Binding energy is the largest source of 
this lowering. The most probable value of the E[^^ dis- 
tribution is made to coincide with i^truc = 235.5 MeV 
with the addition of 45 MeV: E^^^ = -B^oc+45 MeV. This 
correction constant, applied to all z^e candidate events 
at all energies, is irrelevant to the oscillation sensitivity 
reported here-it is merely added so that reconstructed 
neutrino energy and true neutrino energy can be directly 
compared. Of course, the shapes of the signal and back- 
ground i^icc distributions are critical for sensitivity to 
oscillations. The actual relationship between Srec and 
i^truc, including all relevant detector and nuclear effects, 
can eventually be understood with the thousands of mo- 
noenergetic i/^ events and many hundreds of intrinsic Vf. 
events expected during data taking. 

The kinetic energy of the outgoing particles in 
235.5 MeV z/g events is shown in Fig.|4] Assuming 100% 
detection efficiency, the expected rate in terms of E'rcc, 



including smearing effects due to detector resolution and 
nuclear physics, can be seen in Fig. [5] The signal dis- 
tribution's two hump shape comes from the nuclear shell 
structure and the neutron energy levels described by the 
argon-specific spectral function. 



about 10 MeV lower than the monoenergetic signal. Con- 
tributions from decay-in-flight K^ and K^ — >■ w 
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FIG. 4. The reconstructed kinetic energy of the outgoing 
electron and reconstructed total kinetic energy of all outgo- 
ing protons in 235.5 MeV electron neutrino events simulated 
with NuWro. Note that detector resolution and nuclear ef- 
fects are included and events with no outgoing protons (due 
to proton absorption in the nucleus) do not enter the proton 
distribution. 
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FIG. 5. Electron neutrino rate with detector resolution and 
nuclear effects included. The vertical lines designate the sig- 
nal region. 



V. BACKGROUNDS 

The main contribution to the u,, appearance back- 
ground is from intrinsic K'^ — >■ ir^e^Vf^ decay-at-rest. 
The maximum neutrino energy for this decay is 226 MeV, 



are also important and are apparent beyond the three 
body positive kaon decay's endpoint. The relevant v,, 
background contributions to the flux can be seen in 
Fig-H 

Although the uncertainty in kaon production is rather 
large, as exemplified by the ^25% discrepancy between 
the MARS and Geant4 predictions, background measure- 
ments are available in-situ when the experiment turns 
on. Background constraints come from the measured K'^ 
three body decay-at-rest v^ spectrum and its well pre- 
dicted shape as well as measurements of decay-in-flight 
i^e background outside of the signal region. Charged kaon 
induced I'e and monoenergetic v^ events can also be used 
to determine the cross section ratio a^, j On^ for precisely 
transforming measured rate into an oscillation probabil- 
ity. 

Electron anti-neutrino events {peV ^ e'^n), usually 
originating from K^ —> 7r+e~I7g, are a negligible back- 
ground compared to the v,, contribution due to (1) the 
difference in cross section and (2) the minimal V^ flux 
from K^ production and subsequent three body decay. 

A simple single bin, counting experiment is employed 
to discern a potential monoenergetic oscillated Vg sig- 
nal from background. All t^e-like events that enter the 
signal energy region of i^rec G 230 — 255 MeV are con- 
sidered to have E^truc = 235.5 MeV for the purposes 
of the oscillation analysis. Approximately 45% of the 
true i^e signal remains after employing this strict recon- 
structed energy requirement. The range has been chosen 
to provide sensitivity to the LSND allowed region with 
a single bin experiment, in consideration of maximizing 
signal-to-background and having enough events left over 
to discern a signal. Further optimization is possible as 
more knowledge becomes available. The actual experi- 
mental signal region can be refined with the thousands 
of expected non-signal reconstructed Vf^ events in and 
near this energy range as well as ~450 i>e events between 
-Eroc € 100 — 230 MeV. Along with signal region tuning, a 
more sophisticated analysis beyond the simplistic single 
bin, counting experiment can be performed for oscillation 
sensitivity improvement. 

With a 2 kton detector, 2x10^^ POT exposure, and an 
oscillation probability of P{v^ — > i^e) = 0.001, approx- 
imately 16 i^e-like (signal+background) events above a 
background of 3 events are expected. The signal events 
are searched for among an abundant set of distinctive 
monoenergetic i^^ events. The 13 true signal events in 
our example can be compared to approximately 9000 t'^ 
events expected in the same energy window. Obviously, 
the experimental concept described here relies heavily 
on the ability to differentiate the two classes of events. 
All future LArTPC-based electron neutrino appearance 
experiments rely heavily on the ability to differentiate 
charged current muon neutrino events from electron neu- 
trinos ones. The differentiation is largely based on event 
topology with the muon and electron producing a char- 



acteristic track and shower, respectively, along with the 
decay and capture topology unique to the muon. Further- 
more, calorimetric techniques can be used to differenti- 
ate niuon-induced gammas, as in the case of a radiative 
muon capture, from electrons. Although a muon-electron 
misidentification with a LArTPC can occur in rare cir- 
cumstances [21J , it is vanishingly rare for a misidentified 
true 235.5 MeV Vfj_ event to reconstruct as a Ve event in 
the tight energy window required to be considered sig- 
nal {Ej-cc G 230 — 255 MeV) . Including the reconstructed 
energy requirement, the electron- muon misidentification 
rate is assumed to be < 10""* and is deemed a negligible 
background. Backgrounds due to cosmic ray interactions 
and the neutral current interaction v^e^ — )• Vf^e^ are 
considered insignificant as well. 

A 25% systematic uncertainty on the background rate 
inside the signal region is used for this study. This 
number can be considered a conservative estimate and 
will ultimately depend on the understanding of detec- 
tor/nuclear effects, charged/neutral kaon production un- 
certainties, and the in-situ background measurements 
available. Even with a conservative background uncer- 
tainty estimate, the eventual extracted sensitivity is lim- 
ited by statistics. 



E 
^10 



1 : 



10-1 : 




3ct sensitivity 
5a sensitivity 
LSND best fit 
5 LSND 90% CL 



10' 



10"' 



10"' 



10"^ 



FIG. 6. The sensitivity to the LSND allowed region achievable 
with a kaon decay-at-rest source in combination with a large 
LArTPC. 



VI. SENSITIVITY AND DISCUSSION 

The figure of merit for this single bin, counting only 
experimental study is sensitivity to the LSND allowed re- 
gion near Am^ ^ 1 eV^ assuming a two neutrino oscilla- 
tion hypothesis. The two neutrino oscillation probability 
is written as 



Pi'^^. 



= sin^(26l^e)sin^ 1.27 



,Am2 L GeV\ 
eV2 km^ETJ 



(1) 

The constant baseline and neutrino energy can be 
inserted into the equation to form P(f^ — > i^e) = 
sin^(20^e) sin^(0.86Am^). The sensitivity curves, drawn 
in [sin^(26'^e), Am^] space and shown in Fig.Kjj are based 
on calculations of fully frequentist confidence intervals 
using the profile log-likelihood method [33 . "Sensitivity" 
is derived from the median upper limit that would be 
obtained by a set of experiments measuring background 
with no true signal [34 . The uncertainty on the back- 
ground, estimated at 25%, is included as a nuisance pa- 
rameter in the calculations. Baseline smearing effects due 
to detector length and spread in neutrino creation point 
are neglected. 

The sensitivity curves show that a 2x10^^ POT expo- 
sure with Fermilab's 8 GeV booster in combination with 
a 2 kton LArTPC could either confirm or refute the con- 
troversial LSND result at 5a across most of the allowed 
region. In the case that a large sin^(20^e) is observed, 
such an experiment may also be able to provide sensi- 
tivity to the disappearance of i^e [sin {2dee)] and/or Vf^ 
[sin^(26'^^)| through the three and two body i^+ decay- 



at-rest, respectively. Given a localized neutrino source 
and an extended detector, observation of the "oscillation 
wave" as a function of L/E is also possible. 

One of the advantages of a low signal, low background 
experiment is that any event that enters the signal re- 
gion can be carefully examined. Candidate signal inter- 
actions can be studied with the wide range of variables 
available for event characterization in a LArTPC, includ- 
ing a complete three dimensional picture of the neutrino- 
induced charged particle tracks [35i. The search for a 
few i^e candiate events among thousands of similar re- 
constructed energy Vf^ events is made less daunting by 
this fact. The experimental concept is analogous to a 
dark matter search with only a handful of signal events 
expected among a large swath of easily distinguishable 
background events. 

There are a number of relatively uncertain assump- 
tions used in this study that are worth pointing out: 

(1) reconstructed neutrino energy resolution, in consid- 
eration of the nuclear physics affecting the initial in- 
teraction, the poorly understood intra-nuclear process 
leading to the eventual final nucleonic state, and the 
estimate for LArTPC proton reconstruction capability, 

(2) background systematic uncertainty, (3) kaon-induced 
neutrino flux, from decay-at-rest and decay-in-flight, and 
(4) neutrino cross section. As discussed previously, the 
estimated uncertainties associated with each of these as- 
sumptions and, more importantly, their effect on sensi- 
tivity will be significantly reduced when the experiment 
turns on and begins taking data. Note that conserva- 
tive estimates for both kaon production and neutrino 
cross section, each predicted with two different Simula- 



tion packages, have been employed. Taking the larger, 
perhaps less conservative, of each of the two predictions 
would naively increase the expected signal (and back- 
ground) by over 50%. 



VII. CONCLUSION 



Although the detector mass and exposure used for this 
study (2 kton liquid argon detector with 2x10^-^ POT 
from Fermilab's 8 GeV booster) can be considered op- 
timistic, the idea of using the monoenergetic f^ from 
charged kaon decay-at-rest for a Vf, appearance search 
is one that should be considered for current and future 
intense >3 GeV proton sources around the world. 



An experiment based on an intense kaon decay-at-rest 
source and a large LArTPC detector could provide a 
definitive probe of the LSND result and neutrino os- 
cillations near Arn^ ^ 1 eV^. The source would also 
be an excellent opportunity for pion/muon decay-at-rest 
based oscillation searches (e.g. jM]), a coherent neutrino- 
nucleus scattering program (e.g. |371 138] ). and experi- 
ments to measure neutrino cross sections relevant for as- 
trophysics (e.g. [22]), among other possibilities. 
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